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Abstract 



We formulate super symmetric low energy dynamics for BPS dyons in strongly- 
coupled N = 2 Seiberg-Witten theories, and derive wall-crossing formulae 
thereof. For BPS states made up of a heavy core state and n probe (halo) 
dyons around it, we derive a reliable supersymmetric moduli dynamics with 
3n bosonic coordinates and 4n fermionic superpartners. Attractive interac- 
tions are captured via a set of supersymmetric potential terms, whose detail 
depends only on the charges and the special Kahler data of the underlying 
N = 2 theories. The small parameters that control the approximation are not 
electric couplings but the mass ratio between the core and the probe, as well 
as the distance to the marginal stability wall where the central charges of the 
probe and of the core align. Quantizing the dynamics, we construct BPS bound 
states and derive the primitive and the semi-primitive wall-crossing formulae 
from the first principle. We speculate on applications to line operators and 
Darboux coordinates, and also about extension to supergravity setting. 
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1 Introduction 



The wall-crossing in supersymmetric theories refers to the phenomenon where certain 
one-particle BPS states [1, 2] disappear from the spectrum as the vacuum moduli or 
parameters are changed continuously. The naive stability argument of BPS states 
relies on the short-multiplet structure due to partially preserved supersymmetry, but 
this is really applicable only when we consider dynamical processes in a given vac- 
uum. When we change vacuum or parameters, even continuously, the state itself can 
disappear from the spectrum altogether at which point the supermultiplet structure 
of the state becomes a moot issue. 

Although the wall-crossing had an early precursor in the context of supersymmet- 
ric kinks in two-dimensional N = (2, 2) theories [3, 4], it was in the context of N = 2 
supersymmetric theories in four dimensions, such as Seiberg-Witten theory [5, 6] and 
Calabi-Yau compactification of type II string theories, that the phenomenon came un- 
der wide scrutiny. The co-dimension one surface across which a BPS state disappear 
is called the marginal stability wall (MSW), and their presence renders the problem 
of finding BPS spectrum extremely complicated. For the simplest of Seiberg-Witten 
theories, monodromy properties [5] alone can determine the spectrum [7] but this is 
more an exception than a rule. 

Despite such early discoveries, the space-time picture of exactly what happens 
to the state upon the wall-crossing remained unclear until it was uncovered in the 
context of 1/4 BPS dyons in N = 4 super Yang-Mills theory, which preserve four 
supersymmetries just as 1/2 BPS objects of N = 2 theories do. It was found in 
Ref. [8] that such BPS dyons must be, typically, thought of as a loose bound states 
of more than one dyonic centers with mutually non-local charges. The distances 
between such centers are not free but determined by the vacuum moduli u^s 



with (<2yi,p^) being the charges of the A-th center. In particular, some of Rab was 
shown to diverge as a MSW is approached; this happens simply because scalar forces 
and electromagnetic forces do not cancel each other between dyons of mutually non- 
local charges and the equilibrium distance is determined by a detailed balance of 
classical forces: state by state, the wall-crossing has a very mundane and classical 
explanation. 

This finding is immediately applicable to weakly coupled N = 2 theories as well, 
because N = 2 theory BPS solitons can be classically embedded to a iV = 4 theory. 
There will be differences at quantum level because the supermultiplet structures (and 
flavor structures) are different, but the above space-time picture of wall-crossing is 
essentially classical and quite robust. 




(1.1) 
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For both N = 2 and N = A theories, this multi-center nature and the subsequent 
wall-crossing were soon elevated to the semiclassical level [9, 10, 11]. The quantum 
low energy dynamics of magnetic monopoles was derived rigorously from the super 
Yang-Mills theories in question [12], and dyons were constructed as quantum bound 
states of monopoles with certain conjugate momenta turned on [13, 10]. What used 
to be the classical orbit size is now represented by the quantum bound state size, and 
is still determined by the vacuum moduli and charges as in (1.1). The size of the 
bound state is again divergent as a wall of marginal stability is approached, across 
which the state no longer exists as quantum and BPS one-particle state. 

In such supersymmetric low energy dynamics of solitons, precise state counting is 
a simple matter of finding bound state wavefunctions or computing index of ceratin 
Dirac operators on the moduli space. For instance, the bound state of a pair of dyons 
of charge 71 + 72 has been constructed and counted when the total magnetic charge 
is a dual root. The degeneracy on one side of a wall [10] #1 can be written as 

|^(7i + 72)| = 2|(7i,7 2 ) I , (1.2) 

where we introduced ^(7), the second helicity trace for the supermultiplet of charge 

7, as 

n = -itr((2J 3 ) 2 (-l) 2 ^) . (1.3) 

A simplest generalization of this is a chain of dyons with nearest neighbor interactions, 
namely (ta, 7s) 7^ if and only if \A — B\ = 1. Whenever such a state exists as a 
quantum BPS state, the degeneracy takes a simple form [14], 

1^(71+72 +73 +•••)! = jri 2 ^'^ 1 )! • m 

A 

They were shown to exist only when each and every one of {'Ja,1a+i) obeys certain 
inequalities defined by the vacuum moduli, which amounts to being on the "right" 
side of several MSWs, basically one for each interacting pair. The formula clearly 
suggests that such states can be constructed iteratively by attaching one kind of 
dyons at a time, already hinting at a simple universal wall-crossing formula. 

The next breakthrough came from N = 2 supergravity analysis by Denef who also 
found the multi-centered nature of BPS black holes and the subsequent wall-crossing 
in the context of attractor flow solutions [15, 16]. The approach gave a universal and 
explicit constraints for the relative positions of charge centers, say, charge 7^ at x A , 

b+a \xb x a \ 



#1 In this note, we take the convention that Schwinger products take values in Z/2. 
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where Z(^a) is the central charge of 7^ and ( T is the phase factor of the total central 
charge Z T = J2aZ(ia)- In supergravity, this simplifies and supersedes the field 
theory results which we abstractly noted as (1.1). 

The wall-crossing for supergravity black hole solutions is again due to a divergent 
distance between the charge centers, which is dictated by long distance classical 
physics, just as as in the field theory soliton picture of BPS dyons: the sign of the left 
hand side of (1.5) is independent of vacuum moduli while that of the right hand side 
can flip the sign as we change vacuum. Clearly at some point where the right hand 
side approaches zero from the positive side, one of the distances has to diverge, beyond 
which the solution can no longer exist. This is most useful since the sizes of the states 
can be found without detailed construction. However, there is no information on how 
a given charge state is split into what charge centers, unlike the explicit constructions 
of multi-center solution /quantum states in the field theory story. 

Although supergravity solutions themselves were not amenable to explicit and 
precise quantum counting, Denef further went on to conjecture general two-body wall- 
crossing formula that extends the above field theory result to arbitrary (magnetic) 
charges [17]. With spin content taken into account, the formula reads, 

fi(7i + 72) = -(-1) 2 I<™>! 2|( 7l , 72) I n( 7l ) Q( 72 ) , (1.6) 

which was later extended by Denef and Moore to the semi-primitive cases [18], cap- 
tured in a generating function, 

53 !2(Ti + 1172) 9" = JT [ 1 _ ( _i ) t»o.V]- ,t »-»«*»» , (1 . 7 ) 

n=0 k=l 

counting the BPS states of charges 71 + n / ~f 2 in terms of degeneracies of states with 
charges 71 and n7 2 . These spurred much activities toward general solutions to the 
wall-crossing problem, and was integrated recently into more general Kontsevich- 
Soibelman's wall-crossing formalism [19]. 

Despite evidences that support the semi-primitive wall-crossing formulae of Denef- 
Moore (which in turn support Kontsevich-Soibelman formalism), it has been rigor- 
ously tested only in specific cases. The most systematic example of this is the N = 2 
weak coupling analysis that preceded the conjecture, but the limitation of weak cou- 
pling limit casts some shadows on its general usefulness. It would be very useful if 
we can find a similarly systematic method of constructing and counting BPS bound 
states, and apply to diverse BPS objects, such as those dyons that appear in the 
generic strongly coupling region of Seiberg-Witten theory. In this paper, we wish to 
initiate a new framework that can count and construct BPS states, without referring 
either to specific subset of charges or to weak electric coupling, but applicable to a 
large class of iV = 2 theories and BPS states thereof. 
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One common lesson from earlier studies of multi-centered BPS states is that non- 
Abelian completion of the state at charge centers is not essential for understanding 
wall-crossing, since the latter is essentially a long distance phenomenon from the 
spacetime viewpoint. The supergravity solutions are all Abelian while, for solitons, it 
is the long range Coulomb-type interactions that determined the multi-center nature 
of the state. Related is the notion of the "framed" BPS state [20]. The main idea 
there was to treat one or more component dyons as an external object, and the re- 
mainders as dynamical object around such a background. This way, one can treat the 
former as the background, in which the latter moves around and sometime becomes 
supersymmetrically bound to the core state. This split of the state into two parts 
can simplify the state construction and counting substantially. 

Inspired by these ideas, we wish to consider dyons moving around purely Abelian 
dyonic background. In effect, we will split the state in question into the heavy "core 
state" of total charge 7 C and light "halo" or "probe" of charge 7^. For our purpose, 
it is the ratio of the two masses that matters, so this can be for instance achieved 
by approaching a singular point where the probe dyon becomes massless. The low 
energy dynamics of the probe dyon is quite natural thing to do there since, precisely 
at such a singular point, the probe dyon would be the lightest particle among charged 
states. However, there are other circumstances where one part become relatively light 
compared to the other, and our framework will apply. 

Another useful fact is that, as far as wall-crossing behavior goes, we only need 
information near the relevant MSW's, away from which the BPS spectrum is contin- 
uous. This allows another small quantity to play with, by taking vacuum very near 
a marginal stability wall. As we will see later, the distance to the MSW plays a role 
very similar to the weak electric coupling in that it controls the nonrelativistic ap- 
proximation. In the end, we find that the dynamics between the core and the probe 
reduces to massive supersymmetric quantum mechanics with two kinds of potentials. 

These two lead us to a new model of low energy dynamics for dyons in the strongly 
coupled region of N = 2 field theory. Although similar in spirit to the old moduli 
dynamics of solitons, an essential difference here is that the small electric coupling 
constant is no longer needed; this is what allows us to apply the technique to much 
wider class of BPS states than previously possible.* 2 The quantum mechanics has 
four supercharges, as required by the BPS condition, but comes with only 3n bosonic 
coordinates, three for each probe dyon, and 4n fermionic coordinates. Compared 
to the conventional moduli dynamics of weakly coupled regime, we are missing one 
angular collective coordinate for each dyon. This has something to do with the fact 
that we start with dyons, rather than monopoles, as basic building blocks. 

#2 In fact, it should be possible to extend this framework to include gravity and discuss quantum 
bound states of charged BPS black holes. 
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With this new low energy dynamics in place, we can compute how many BPS 
bound states of the core and the probe dyons can form, and under what condition. 
At the end of day we derive, via a first-principle computation, the semi-primitive 
wall-crossing formula with 71 = 7 C and 72 = 7^. In this note there is in fact no 
restriction on 7,, as far as such a state actually exist and all of its component dyon 
centers can be made heavy. Thus we in effect are computing £l(j c + 717^) with the 
only restriction that the dyon 7^ is primitive and become massless somewhere in the 
vacuum moduli space. We wish to emphasize that, alternatively, we may think of the 
theory as a setup for finding framed BPS state with line operator of charge 7 C and 
halos 7 h [20]. 

This paper is organized as follows. In Section 2, we write down the long-distance 
Abelian form of the core state in terms of the central charge function, while the probe 
dyons are treated as quantized solitons in that background. As a result we find a 
bosonic low energy Lagrangian of the probe dyons purely in terms of quantities that 
can be constructed out of the central charge functions. This reproduces some of gen- 
eral results, such as distances between two charge-centers, obtained from supergravity 
attractor flow analysis, even though we are dealing with field theory states. 

Section 3 discusses how one can construct a Af = 4 supersymmetric Lagrangian 
with 3n bosonic coordinates and 4n fermionic coordinates, by extending previous 
studies by Coles and Papadopoulos [21] and also by [22]. These previous works con- 
structed massless supersymmetric theories of similar kind, which is, however, missing 
the crucial elements of potentials. Without the latter, the bound states we are in- 
terested in cannot form at all. We construct in particular massive theories in which 
degrees of freedom are cataloged by SO(A) R = SU (2) L x SU (2) R algebra with bosons 
in (3, 1) (thus, the first SU{2) L also serves as a rotation group) and fermions in (2, 2) 
representations. The four supercharges are also in (2, 2). The Lagrangian has SU (2) R 
symmetry manifest while SU (2) l can be explicitly broken by the background. 

Section 4 shows how the general discussion of section 3 makes contact with the 
probe dyon dynamics of section 2 under the assumption the vacuum moduli of the 
underlying N = 2 theory is very near the MSW. The latter assumption controls 
the energy scale of the potential energy, and allows a nonrelativistic approximation 
possible. We then quantize the resulting dynamics and derive the bound states for 
7c + Ihi and again shows how the bound state size diverges as one approach MSW 
and how the bound state is impossible on the other side of MSW. 

Section 5 elevates this to a primitive wall-crossing formula, and extends further to 
the cases of 7 C + n / ~ff l by invoking spin-statistics theorem. This derives, in particular, 
the semi-primitive wall-crossing formula from a first principle computation. 

We then conclude in Section 6 with summary and other comments especially on 
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how one can make use of this formalism to compute the line operator expectation 
values and how one can extend the formalism to the supergravity setting. Some 
computational details are summarized in Appendices. 

2 Classical Dynamics of Probe Dyons 

In this section, we construct the semiclassical form of the core state, entirely in terms 
of the central charge function, and describe energetics and dynamics of a probe dyon 
in the core state background. This leads us to a bosonic Lagrangian of the probe 
dyon, which will be supersymmetrized and quantized in the later section. 

Although the exercise here applies to any core state one can imagine, as long as 
there solve the relevant semi-classical BPS equation of the effective Abelian theory, 
we are eventually interested in core states that actually exist as quantum BPS states. 
It is known that the former does not always imply the latter [23, 10]. Alternatively, 
for the framed BPS states, the core state should correspond to a supersymmetric 
line operator. Either way, we are interested in case where the supersymmetric lift 
of this probe bosonic dynamics would make sense in the context of the underlying 
four-dimensional theory. 

2.1 Semiclassical Core State 

We start by recalling semiclassical properties of N = 2 dyons when expressed in terms 
of the low energy theory of Seiberg and Witten. Traditionally the smooth solitons 
are possible only when we include the entire non-Abelian origin, but this is practical 
only in the weakly coupled limit. 

To avoid such restrictions, a more convenient starting point is to write the BPS 
equation in the Abelian low energy description of Seiberg and Witten. This approach 
was investigated previously [24, 25] with emphasis on split flow picture of the classical 
soliton and gave an interesting parallel to the string web picture [26] of N = 4 1/4 
BPS dyons. These solutions are invariably singular at the charge centers, since there 
is no non-Abelain mechanism to stop the Coulomb-like divergence at origin, which 
was controlled ad hoc by introducing UV cutoffs. 

For our purpose, however, this divergence is of little consequence, essentially be- 
cause we will be using this solution as background. As long as we can ascertain 
existence of quantum state of such a charge and as long as we put correct boundary 
condition at such singular points, forcing the probe dyon wavefunction to vanish there 
fast enough, there would be no physical problem associated with it. It is entirely anal- 
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ogous to the Hydrogen atom problem of undergraduate quantum mechanics, where 
finite and trustworthy bound states are obtained even though the Hamiltonian is 
naively singular at origin. 

Using the SUSY transformation rule for gaugino along the particular direction 
parameterized by a phase factor (, one can obtain the BPS equations 

f. _ i^Vfa = , (2.1) 

where % labels the unbroken U(l) gauge groups, and J 7 denotes the complexified field 
strength 3- vector B + iE. See appendix A for details. There is also an electric version 
of this equation 

r D ~tC l V ( j ) \ y = 0, (2.2) 

—* . . . —* 

with T l D = r %3 Tj and 

r y = ^r^sw(0) , (2.3) 

where F sw (0) is the Seiberg-Witten prepotential of the given theory. Since it is ReJ^ 
that enters the Gauss constraint, the field strengths are such that [24] 



Re / Ti = AnP\ Re / F D = -4vrQ\ (2.4) 

with the total magnetic charges P l and the total electric charges Q l 

In particular imagine a semiclassical core state of charges 7 C = (P\ Q l ) = 7 c ,a, 
with 7a = (P\, Q\)i distributed into several dyonic cores at x A , and the field strength 
takes the following asymptotic forms, 




(2.5) 

One can show that ( can be identified as the phase factor of central charge Z COTC of 
this core state* 3 

Zcorc = |Z corc |C = QVi(oo) + Pi^(oo) . (2.6) 
This semiclassical description is, strictly speaking, valid away from x = a^'s. 



#3 Scc Appendix A. 
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Note that the positions, x^s, of the centers would be restricted by an analog of 
(1.5). Precise positions of these centers is, however, immaterial for counting BPS 
bound states, as long as the relevant core state actually exists as quantum and BPS 
bound state. This happens because one ends up computing supersymmetric indices, 
which are robust under small deformations of the supercharges. More important is 
how the core electromagnetic charge is distributed into such centers. See section 4 
for related discussions. 



2.2 Probe Dyons and Electromagnetic Forces 



Let us now introduce a probe particle of charge 7^ = {pi,qi), in a background cre- 
ated by such a core state. It will be considered as a probe particle in the external 
electromagnetic field by the massive core state. Using the equations (2.1,2.2), one 
obtains 



q-f + p-f D = iC 1 VZ h 



(2.7) 



where Zh = q- (fi+p- 4>d is now understood as position-dependent. We introduced the 
notation Z to emphasize that this quantity is position-dependent. The usual central 
charge Z is related to it as Z = Z(oo). 

The real and imaginary part of the relation will give us hints how to construct the 
low-energy Lagrangian of probe dyon in the background of core particle. The real 
part can be succinctly written as 



VVbouiomb = -VRe 



(2.8) 



where 



Vbouiomb = Re ( r ) i j ^2 



P r A 



X 



xa\ 



+ (Re(r))- 1 E 



($ + Imtfijpi) (Qj,A + Im(r)ij-Pj) 



\x — Xa\ 



(2.9) 



is nothing but the Coulomb potential energy felt by the probe dyon due to the core 
state. The real part of this equation is even simpler 



E 



Qa-P- Pa- q 



-Vim 



\x 



XA 



C l z h 



or equivalently 



(lc,A, lh) 



-Vim 



\x 



XA 



C l z h 



(2.10) 



(2.11) 
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which, as we will presently see, encodes the Lorentz force on the probe dyon. #4 

We first discuss the invariant expression of the minimal coupling under the Montonen- 

— * — * 

Olive duality. Recall that, from the BPS equations, one can conclude that (J 7 , Fd) 
transform under the duality transformation as vector representation like (0, 4>d)- For 

— * — * 

example, let us consider the S-duality transformation of (J 7 , Fd) 

B ->■ -Im(r)S + Re(r)B , E ->■ lm(r)B + Re(r)E . (2.12) 
Then, one can easily show that, under the S-duality transformation, 

q -> p , v ->■ -q , (2-13) 

where we used, for the last transformation, the fact that r — > — r _1 . 

When the probe dyon moves (slowly) under the electromagnetic field of core par- 
ticle, the minimal coupling terms therefore become [27, 28] 

£ int = q*v ■ A* + tfv ■ A % + q l A\ + p% , (2.14) 
which is the duality invariant expression. Here A^ and A^ are defined as 



R e jr = V x A , ReTo = V x A , 
Imf = V • A , InuF D = V • i ■ 



(2.15) 



Using the BPS equation (2.7), the interaction terms can be managed into a rather 
simpler form 



C 



mt 



-v-W + Re 



C l Z h {x) 



Re 



C _1 ^(oo) 



(2.16) 



where the vector w satisfies the relation below 



V x W = Vim 



(2.17) 



Note that, in (2.16), Re[C 1 Z h (oo)} = Re[C 1 Z h ] represents the lowest possible energy 
the probe dyon can attain. 



#4 The normalization of charges and the sign convention for Schwinger product differs from that 
of Rcfs. [15, 18] 

(7:7') = ^^'^Dcncf-Moorc 
Related is the fact that Olir £ is CDenef-Moorc- 
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2.3 Massive Moduli Dynamics of Probe Dyons 

Finally we come to the effect of the long range scalar field on the dyon. The low- 
energy Lagrangian of probe dyon 7^ moving in the background of core particle 7 C can 
take the following form 

£ bosonic = £k . n + £ . nt ^ (218) 

where the kinetic term must be [27, 28] 

Adn = -\Z h (x)\VT^~ -\Z h (x)\ + l -\Z h {x)\v 2 + 0{v A ) (2.19) 

with Z^ix) = q ■ <p + p ■ </)£>, since |^(x)| is the effective inertia of the probe dyon. 
Adding all these together, we find the classical Lagrangian, 



: hosonic = Uz h (x)\v 2 -\z h (x) 



Re {C l Z h (x)) - ReJr^ftM] -v-W (2.20) 



with V x W = Vim (C^Zhix)). 

This Lagrangian has the classical ground state at x = x* where [^(rc*)] = 
Re[C _1 2/i(a;*)], with the ground state energy Re[C _1 >2/ l (oo)]. We wish to elevate 
this, later, to M = 4 quantum mechanics, so it is more convenient to separate out 
the ground state energy. Thus, our starting point is the bosonic Lagrangian, 



/'bosonic /'bosonic , p 

^moduli _ i- rte 



C L Z h (oo) , (2.21) 



so that supersymmetric bound states would have zero energy. This also reproduces 
an analog of Denef 's formula [15] for the probe dyons since, 

J2i^ i: ^ = HC 1 Z h (oc)} . (2.22) 

\Xa 

from Eq. (2.11) and Im^" 1 Z h (x *)] = 0. This is the same as (1.5) once we realize 
that total central charge Zt = Z c + Zh is dominated by Z c since Zh/Z c is very small; 
(t is approximately equal to (. 



2.4 Fermionic Partners 

We have derived a classical (thus purely bosonic) Lagrangian that describe the dy- 
namics of a probe dyon in the background of the core state, with 3 bosonic collec- 
tive coordinates per each probe dyon. Without much effort, we can further deduce 
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that each probe dyon will also come with 4 fermionic degrees of freedom, giving An 
fermionic variables as opposed to 3n bosonic variables. 

The simplest way to see those four fermionic variables is to recall that a BPS 
particle, of a given charge, in N = 2 theory are at least in the half-hypermultiplet, 
with spin content 

[1/2] ©2[0]. (2.23) 

This spin content can be generated only if the dyon comes with a pair of complex 
fermionic degrees of freedom in a spin 1/2 multiplet, which translates to four real 
fermionic coordinates. They are, when we consider the dyon in isolation, also Gold- 
stino modes coming from the four supercharges broken by the BPS state. More 
generally, the probe dyon could be in a BPS multiplet of type, 

[s]®([l/2]©2[0]) , (2.24) 

with an angular momentum multiplet [s] of spin s, in which case [s] typically arises 
because the probe dyon is itself a composite or has, otherwise, some internal light 
degrees of freedom. What matters for us is that we still have the same four fermionic 
collective coordinates whose coupling to the bosonic ones are tightly constrained by 
the M = 4 supersymmetries. 

When we consider the special limit of solitonic dyons in weakly coupled theo- 
ries, this mismatch between the bosonic and the fermionic degrees of freedom can 
be understood easily [29, 30, 37]. Solitonic dyons arise there from excitation of a 
monopole soliton with particular U(l) momenta turned on [31]. While the initial 
monopole soliton comes with four bosonic and four fermionic collective coordinates, 
one angular bosonic coordinate is traded away in favor of its conjugate momentum 
(which is physically the electric charge). This procedure, however, leaves the four 
fermionic coordinates intact. It has to be so, since the dyon is still BPS and the 
necessary half-hypermultiplet structure would be generated using all four of these 
fermionic degrees of freedom. Nor does this reduce the Af = 4 supersymmetry of 
the remaining dynamics, although their embedding into the underlying field theory 
is rotated in response to the new electromagnetic charges. 

3 Massive M = 4 Mechanics onto Moduli Space 

An odd fact, when we consider a supersymmetric lift of the above Lagrangian for 
probe dyons, is that the low energy dynamics involves 3 bosonic collective coordinates 
for each probe dyon, yet, there should be 4 fermionic counterparts. Supersymmetry 
with mismatching bosonic and fermionic degrees of freedom is in principle possible 
for quantum mechanics because there is no notion of spin, but still construction of 
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such theories, especially with extended supersymmetry, was not widely studied. The 
only known example is certain (massless) class of supersymmetric nonlinear sigma 
models by Coles et. al. [21], which were later specialized in the context of extremely 
charged black holes of the same charges [22]. Neither of these studies considered 
massive versions, as needed here, however. 

Similar situation existed a dozen years ago when low energy dynamics of solitonic 
monopoles were studied for weakly coupled iV = 2, 4 Yang-Mills theories. The con- 
ventional massless moduli dynamics [32, 33, 34] with 4n- dimensional target manifolds 
without potential were found to be inadequate for dyons in generic Coulombic vac- 
uum when the rank of the gauge group is two or larger [8]. The problem was the lack 
of potential terms in this older formulation. The low energy dynamics of monopoles 
had to be reformulated so that both the potentials and Af = 4 supersymmetry are 
manifest. Later, such massive Af = 4 quantum mechanics were found, simply by 
twisting supercharged by triholomorphic Killing vector fields on the moduli space 
[9, 11, 10, 12]. * 5 This lead to a whole machinery whereby dyon spectra in the weakly 
coupled limit of N = 2,4 Yang-Mills theories were constructed explicitly [14]. See 
Ref. [37] for a broad overview of this development. 

In this section, we wish to investigate how the new kind of classical low energy 
dynamics of section 2 can be also elevated to one with Af = 4 supersymmetry. We 
will find that massive Af = 4 supersymmetric mechanics with mismatching bosonic 
and fermionic degrees of freedom is possible and will, specifically, build a massive 
(i.e. with potential) supersymmetric Lagrangian with 3 bosonic coordinates and 
4 fermionic coordinates. This restriction to the lowest possible target dimension 
simplifies the construction greatly, in part because the target manifold turned out to 
be conformally flat R 3 , and yet still good enough for deriving semi-primitive wall- 
crossing formula.* 6 

3.1 Toy Model: Flat R 3 Target 

As a toy model, let us pretend that the bosonic moduli space is flat R 3 and see how 
scalar and vector potentials on R 3 can be incorporated into the quantum mechanics 
in a manner consistent with four supercharges. 

Af = 1 supersymmetry is easy to incorporate. We start with the usual transfor- 

#5 Some related mathematical structures were first studied in Rcfs. [35] while its potential connec- 
tion to dyons was previously hinted by Ref. [36]. 

# 6 For generalization that can address many probe dyons with non-negligible mutual interactions, 
we need to consider higher dimensional target manifolds, which is left for a future work. 
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mat ion rule 

§ x a = -ie^ a , 5i> a = ex a , (3.1) 
under which the following free Lagrangian that is invariant 

= l -± a x a + % -^ a i) a . (3.2) 

Since we are dealing with quantum mechanics, rather than a field theory, we can add 
any number of fermions, as long as we let them be invariant under the above super- 
symmetry transformation. As we will see shortly, however, extended supersymmetry 
would not leave this extra fermion intact. 

For our purpose, one extra fermion A is needed for each triplet of (x a ,ip a ), so we 
may start with 

£ (0) = \x a x a + + ^AA , (3.3) 

where 

5X = . (3.4) 

Incorporation of an external gauge field w on R 3 is equally easy. Adding a minimal 
coupling —x a w a to the Lagrangian and noting the supersymmetry transformation 
property, 

5( — w a x a ) = + ieil} a x b (d a w b — d b w a ) + total derivative 

= + id a w b {ei> a x b - ei> b x a ) , (3.5) 

one finds a canceling term of type 

5( + id a w b i> a i> b ) = +id a w b {x a ei) b - x b eip a ) . (3.6) 

In summary, the following Lagrangian has Af = 1 supersymmetry 

C = l -x a x a + ^ifj a tp a + ^AA - w a x a + id a w b ifj a ip b . (3.7) 

In order to introduce the bosonic potential to the above model, we modify the trans- 
formation rule of the auxiliary fermion A as 

5x a = -ie^ a , 5^ a = ex a , 5X = eK , (3.8) 
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upon which we find 




(3.9) 



The canceling term for this is 



5( + id a Kij a \) =ie\x a d a K - ie^j a Kd a K 
=ie\x a d a K + 5x a Kd a K , 



(3.10) 



while one must add one more to close the transformation algebra, 




(3.11) 



In summary, the Lagrangian 




- x a x a + l -rr+ l ^-wax' 




(3.12) 



has Af = 1 supersymmetry, for any K and w. 

We eventually wish to formulate dyon dynamics with J\f = 4 supersymmetries. For 
conventional supersymmetric quantum mechanics, this requires the target manifold 
to be 4n dimensional and hyperKahler, which is clearly inappropriate for our 3n 
dimensional target. Nevertheless, the BPS nature of the dyons and existence of BPS 
bound states implies that there should exist such an M = 4 lift. 

To find the relevant supersymmetries and the subsequent restrictions on the po- 
tentials, note that, since the number of bosons and the number of fermions mismatch 
by 3 to 4, we can organize the degrees of freedom using 5*0(4) = SU(2) L x SU(2) R 
algebra. Let the bosons, x a ,a = 1,2,3, transform as (3, 1) representation while the 
fermions, (ip a , A), are naturally in (2, 2) and better denoted as ip m , m — 1, 2, 3, 4 with 
ip A = A. Thus, a,b,. . . are the vector indices of SU{2)l while m,n,... are vector 
indices of SO(A). The Af = 4 supersymmetries are then naturally in (2,2) under 
this SO (A), since it should relate bosons to fermions. Thus the four supersymmetry 
transformation parameters will be denoted by e m . 

A useful method of relating SU(2) L objects to SO(A) object is to employ 't Hooft's 
self-dual symbol Based on previous experience of embedding SO(3) ~ SU(2)i 
into SO (4), such as in Yang-Mills instanton construction, one can guess the following 
M = 4 SUSY transformation rules 




a n-a , 7^ 



(3.13) 



(3.14) 
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which, for e 4 = e, matches (3.8). 

This suggests that the Lagrangian (3.12) can be extended to admit M = 4 super- 
symmetries, if we can organize the fermion bilinears in terms of 77 symbol as 

1 3 i 4 i 1 

C = -^x a x a + ^V™ + 2^WT - w a x a ~ ^K 2 , (3.15) 

a=l m=l 

which matches (3.12) if we impose 

e abc d a K = d b w c - d c w b . (3.16) 

One can indeed show that the above Lagrangian is invariant under the M = 4 SUSY 
transformation rules (3.13). This Lagrangian is manifestly invariant under SU(2)r. 
The SU(2)l invariance is broken only to the extent that K breaks the rotational 
invariance. If K is spherically symmetric, for instance, the full 5*0(4) symmetry 
would be restored. 

Let us discuss the closure of the M = 4 algebra. For bosonic variables, one can 
show 

5 c 5 e x a = - tv * mnV b pn e m Cx b + t V a mn e m CK , (3.17) 

which implies 

(5 c 8 e - 8 e 5 c )x a = -2ie m ( m x a . (3.18) 

Here we used the following identity r/mnVpn = <5 a6< W + ^Vmp- Let us now in turn 
consider the case of fermionic variables. 

(5 ( 5 e - 5 £ 5 c )^ m = - 2ie n C^m + i(e n ( m + e m ( n )r + ir] a pq d a K {e\ m + e m C)^ , 

= -2ie n C^m , (3-19) 

where for the last equality we used the equation of motion of ip m 

r + v a qn daKtp n = . (3.20) 

We therefore conclude that the Af = 4 SUSY algebra is given by 

{Qm, Qn} = 25 mn H , (3-21) 

with the Hamiltonian H. 
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3.2 Toy Model with M = 1 Superfields 



We can write the above Lagrangian by introducing M = 1 superspace with an anti- 
commuting coordinate 9. Following the notation in Ref. [22], we define the bosonic 
and the fermionic superfields as 

= x a - ieip a 7 A = i\ + i9b . (3.22) 

The supersymmetry generator and the supercovariant derivatives are then, 

Q = d e + iOdt, D = d e - idd t . (3.23) 

Our toy model based on flat R 3 , with scalar and vector potentials, can be written in 
a superspace form as 

C = J d6 (^-D<$> a d t <$> a -^ADA + iK(<f>)A-iw(<S>) a D<f> a ^J . (3.24) 

Although only M = 1 supersymmetry is manifest, we saw that Af = 4 supersymmetry 
will emerge if the condition *dK = dw is imposed. This form of the Lagrangian is 
useful because it could be generalized to the curved moduli space almost immediately. 



3.3 Massive M = 4 Theory onto Conformally Flat R 3 

Recall that, for a single probe dyon, there are three quantities that appears in the 
bosonic moduli dynamics. The scalar and the vector potentials, as we already in- 
corporated into M = 4 toy model above, and most crucially, the position-dependent 
mass term \Z h \ for the coordinates x a . Thus, in addition to the above interaction 
terms, we wish to replace R 3 by a conformal flat R 3 whose metric is 

g a b = f5ab , (3.25) 

with / later to be identified with \Zh\- In fact, as can be inferred from the massless 
version in Refs. [21, 22], M = 4 supersymmetry restricts the three-dimensional metric 
to be conformally flat. 

We defer detailed construction to appendix B, and simply state here that the 
desired Lagrangian, now with potentials, has the superspace form 

c = j de (^-f^)D^ a d t ^ a - \f(®)hDA 

+ie at A/($)D$ k S$ c A + *($)A-iW($) a Dfj , (3.26) 
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with the condition 

d a IC = e abc d b W c (3.27) 
imposed. In terms of component fields, this equals 

C = ^f(x a x a + i^ m V t ^ m ) 

' {2d 2 f-r\df) 2 )e mnpq ^ m rW 



4-4! 



-^/C 2 - W a x a + \j xl2 d a {r Xl2 K)rLJ> m r , (3-28) 

where the covariant derivative for fermions is given by 

= i, m + ^e^f-'dj^r . (3.29) 

As in the flat case, the degrees of freedom and the supercharges are cataloged by 
50(4) = SU{2) L x SU(2) R algebra, and the Lagrangian is manifestly invariant under 
SU(2) R . The £77(2) L keeps track of how / and JC (and thus W also) transform under 
spatial rotations, and become a symmetry whenever these quantities are rotationally 
invariant. 

This £0(4) structure and SU(2) R symmetry tells us an extended M = 4 super- 
symmetry exists, as in the flat R 3 example. It is not difficult to see that 

6 e x a = iv a mn e m r, 

= rf mn e n x a + e m b, (3.30) 

with four Grassman parameters e m leaves the Lagrangian invariant. The only differ- 
ence from the flat case, (3.17), is that K is replaced by its generalized form, namely 
on-shell value of the M = 1 auxiliary field b, 

b = 7 ( K + ^pM^ q ) ■ ( 3 - 31 ) 

The superalgebra remains the same as the flat case, 

{Qm,Qn} = 15 mn H , (3.32) 

where we denoted the four supercharges by Q m as before and the Hamiltonian by H. 
For completeness, we also record the classical form of the Hamiltonian, 

H dassical = ^V + ^{2d 2 f-r\df) 2 )e mnpq i, m rri> q 

+ ^K 2 - ^.(r^Cff , (3-33) 
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with the covariantized momenta 

r = Pa + W a - l -e ahc d h frf m J, m r . (3.34) 

The quantum Hamiltonian differs from this by normal ordering issue, and can also 
be found in appendix B. 

4 Quantum BPS States near WMS 

4.1 J\f = 4 Low Energy Dynamics of Dyons near MSW 

Let us stop here and ask under what circumstances we actually expect to see a sensible 
low energy dynamics of dyons to appear. The old setting based on dyons as quantum 
bound states of excited magnetic solitons was possible by resorting to weakly coupled 
regime. There, the basic requirements was that the energy due to electric charges 
and also due to motion of the solitons are of higher order. Thus, the reduction to 
quantum mechanics is controlled two small parameters; typical speed of the magnetic 
soliton and the electric coupling constant [32] . 

Here, however, we are here dealing with dyons of generic charges at generic cou- 
pling, and must find different criteria to justify reduction to low energy quantum 
mechanics. Note that the weak coupling requirement and the small speed require- 
ment of old moduli dynamics is in fact interrelated. That happens was that the 
moduli dynamics of N = 2 and N = 4 monopoles usually acquire a bosonic potential 
of order e 2 , so for typical states the small electric coupling is necessary to ensure 
small velocities. 

In the present low energy dynamics of probe dyons around a core state, the size 
of the potential is instead controlled by how far are the phases of central charges 
of core and halo particles are aligned. Thus, by staying very near MSW, we have 
a good control over the potentials. Furthermore, the massgap between this sector 
and the rest is also substantial, and controls possible interference from other charged 
particles.* 7 So it is the proximity to the MSW and also the mass ratio of the two 
parts that now control the reduction to the low energy quantum mechanics. 

With this mind, we compare (2.21) against the supersymmetric Lagrangian (3.28). 

#7 The latter is easiest to see when the small mass ratio is achieved by being near a singular point 
of the vacuum moduli space. The relevant coupling that governs the interaction of the held theory 
would be a dualized coupling which becomes small as the singular point is approached. 
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One can see the supersymmetric uplift may work only if 



f=\Z h \ , 



^-K 2 = \z h \-Re[c 1 z h ], 



V x W = V (imfC 1 ^]) . (4.1) 



but the requisite Af = 4 relationship between /C and W, *d/C = </W, is not yet 
apparent. Thankfully, this condition is satisfied precisely when the criteria for the 
low energy approximation are met, as we described above. 

To see the latter, write C^ x Z h = | Z h \ e t(3 . Near the wall of marginal stability, the 
angle (3 at spatial infinity is very small whereas its value at classical vacuum is 0. 
Recall that the bound states we wish to find and count are all peaked at the classical 
vacuum manifold. This allows us to expand relevant quantities in small (3. As we 
move closer to charge centers, f^'s, (3 can grow again but the precisely form of the 
background at such charge centers are not to be trusted and also happily immaterial 
for our purpose of finding BPS bound states. Therefore, we take the value of f3 to be 
small everywhere and find 



K 2 = 2| Z h \ 2 (l - cos/3) ~ | Z h \ 2 f3 2 ~ | Z h \ 2 (sin(3) 2 = (imfC 1 ^]) 2 . (4.2) 



Thus, for all practical purpose, we may identify /C = Imf^ -1 ^] and the Af = 4 
requirement (3.27) is obeyed automatically. This completes the derivation of Af = 4 
moduli dynamics in (3.28) of a probe dyon in a given core state background. 

The function /C can be generally written, from (2.11), as 



with 7 Ci ^4 centers of the core states at xa and also /Co = Im[^ _1 ^(oo)]. Details of 
/ = \Z h \ won't matter much for the purpose of constructing bound states, it turns 
out, as long as we keep track of its singular behaviors at charge centers. 

Before we start the detailed analysis, let us note again that the semiclassical core 
state here is not really a good representation very near its charge center (s), where the 
non-Abelian nature of the states becomes relevant. Naturally, the low energy dynam- 
ics of probe dyons is plagued by the same issue. However, this hardly matters near 
MSW because the bound state (if any) would be very large and determined entirely 
by Abelian part of the low energy field theory: Whatever singularity at Coulombic 
centers cannot alter such wavefunction significantly, as long as we impose the bound- 
ary condition at centers intelligently enough. This should become more obvious when 
we discuss actual bound state wavefunctions in section 4.3. For this reason, we may 
as well take the above form of /C, /, etc literally, and consider supersymmetric bound 
states thereof, with some care given to the boundary condition of the wavefunctions 
at centers Xa- 




(4.3) 
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4.2 Quantization and Supercharges 

Let us start with the canonical commutators. The conjugate momenta of bosons are 
denoted as p a , 

[ Pa ,x b ] = -t5 b a , (4.4) 

whereas the normalized fermions, tfj m = f 1 ^ 2 ^ 171 , are more convenient for writing out 
the remaining canonical commutators, 

^} = 5 mn , [ Pa , 4> m ] = = [x a , V>] . (4.5) 

With this we can now write the four supercharges as 

Qm = ~V a mn r(Pa + W a ) + ^L/'^Jf + ^./^"W^ + JC^ . (4.6) 

For the proof that these are right supercharges, see appendix B. In particular, the 
supercharge associated with e 4 is 

Q = Q A = r (p a + W a ) - \r l d a fr + \d a fe ahc ^^ + A/C . (4.7) 

Since the superalgebra implies {Q m , Q n } = ^mnH, the ground state of the system 
can be found by demanding that it be annihilated by Q4. 

4.3 BPS Bound States and Marginal Stability 

The canonical commutator of the fermions 

is a Clifford algebra which can be represented by Dirac matrices, 

v / 2^ = 7 a =(° a o"), v/2A = v/2^ = 7 4 =(_°, j) , (4.9) 

and wavefunctions can be regarded as 4-component spinors on R 3 . Also useful is the 
chirality operator 

r = yyyy = ( J \ ) • (4.10) 

Under the above representation, one of supercharge Q4 now has a simple form, 

VV Q4 = l^Pa + W a ) ~ '-(dahgfh 11 ^ + IC^ , (4.11) 
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or more explicitly, 



PTfO - I ° a-(p + W)+ilC \ 

VZJ \a-( P + W)-ia-d(\ogf 1 / 2 )-ilC )• 

We wish to find supersymmetric ground states, = 0. Since "H^ = then, such 
states would actually preserve all four supercharges. Such states are then automati- 
cally BPS with respect to the N = 2 field theory with the energy Re[C _1 ^(oo)], as 
can be seen from (2.21), not counting the core state energy. 

Write the four- component wavefunction \l> as 

* = ( r 7 ) , (4.13, 

upon which two component wavefunctions U and V obey 

{<r-{p + W)-iK)U = 0, 

(cr • (p + W) + iK) V = 0, (4.14) 

With the supersymmetry condition dJC = *dW, it is easy to see that the first equation 
cannot have a normalizable solution while the second may. Denoting the respective 
operators as T>±, 

V T V± = (p + Wf + IC 2 + a a {d a K ± d a K) , (4.15) 

which shows that V + V_ is a positive definite operator while V_V + is not. Only the 
latter can have zero modes. Thus, we arrived at the conclusion that the counting 
of BPS bound states between the core dyon and the probe dyon becomes that of 
counting normalizable two-component zero modes V of the operator T> + , with the 
final form of the BPS bound state 

*=({0' (4 - i6) 



with V+V = 0. 



It is illuminating to solve this equation for the particular case of spherically sym- 
metry core state. The vector potential would be that of a Dirac monopole, so we 
denote 

W = -gA Dirac , g = -(7c, 7h), A Dirac = -cos9d(fr , (4.17) 
from which follows the scalar potential 

K = /C + - . (4.18) 

r 
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In this case SU(2) L also becomes a symmetry, allowing an explicit solution to the 
bound state problem. The number g is half-integer quantized, as dictated by the 
Dirac quantization of this quantum mechanics, and also from its field theory origin 
as the Schwinger product of the two quantized charge vectors. The angular and spin 
part of the wavefunction is classified by spinorial monopole spherical harmonic tensor. 
The lowest possible angular momentum would be than j = \g\ — 1/2, since the charge 
interacting with such a Dirac monopole, W, is endowed with a well-known angular 
momentum —gf. Tensoring with the intrinsic spin 1/2, the minimum possible value 
j — \g\ — 1/2 follows. 

Denoting the corresponding the lowest-lying two-component angular momentum 
eigen-states r)j=\ g \-i/2,m of 5*0(3) ~ SU(2)l we rely on Kazama et.al. [39] for reduc- 
tion of the above to the radial equation, 

V = M0»fe=|*l-i/2,m, (- i fg~\ X Tr + \ + MO = ■ (4.19) 

Integrating the latter equation, we find 

h(r) = ^exp(^- J /C(r)^) =Orl^^ 1 exp(-[sgn(( 7c ,7 h ))-/Co]-r) , (4.20) 

with the normalization constant C. Note that this gives a normalizable ground state 
if and only if the half-integer-quantized (7,, 7^) is not zero and has the same sign 
as /C = Im[C _1 ^(oo)]. #8 The latter condition is also reflected on the fact that the 
probability density of this wavefunction is peaked at radial size 

(7c 7.) (421) 



Mc-^oo)] ' 

which, for a single-center core state, exactly mirrors the classical orbit radius in (2.22). 

The sign of (7 C , 7^) is determined by the charges of the core state and the probe 
state, and does not change as we move along the vacuum moduli space. However, the 
/Co = Im[C _1 Z/ l (oo)] does change its sign across the marginal stability wall between 
the core state and the probe state. Classically, this happens because glC < would 
make the potential repulsive. The upshot is that the BPS bound states of one side 
where (7 C , 7/ l )/Im[^ _1 2/ l (oo)] > disappear as we move to the other side where 
(7c 7/i)/ im [C~ 1 ^(°°)] < 0, as was originally found in the supergravity setting. 

With this exercise, we learned a few things: 



# 8 L 2 normalizability requirement from r = region is satisfied as long as \(g c ,Jh)\ is not zero, so 
does not impose additional restriction. 
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• Normalizable bound state between the core state and the probe state is realized 
only when the Schwinger product of the two charge is nonzero. 

• Normalizable bound state between the core state and the probe state is realized 
only when the Schwinger product of the two charge is of the same sign relative 
to the value of Im[^ -1 ^] at spatial infinity. 

• When such normalizable states exist, the degeneracy is 2j + 1 = 2|(7 C , 7 h )|. 

Much of the above statements are properties of a Dirac operator with V± as the chiral 
and the anti-chiral parts; there must be an index theorem associated with them. 

In fact, the structure of the operators are essentially that of an electrically charged 
fermionic field around the magnetic monopole, except that we do not see the non- 
Abelian structure that regulate the short- distance behavior of the core state. Similar 
issues in the context of quantization in the backgrounds of non-Abelian monopoles 
vs. Dirac monopoles (or more precisely Wu-Yang monopoles [40]) have been studied 
in depth decades ago, where it was found that with proper boundary condition at 
origins of the latter, behaviors of the two are essentially the same [41]. The boundary 
condition is constrained by the requirement that the Dirac operator constructed out of 
T>± should be Hermitian, which is known in the literature as the self-adjoint extension. 

This is related to the fact that, even though the two potentials of the quantum 
mechanics are singular at origin, the wavefunctions found are regular everywhere and 
in particular suppressed strongly at origin. If we attempted to solve for T)JA = 0, the 
radial eigen-function of U would have the behavior r - ^ 70 ' 7 ^' -1 at origin and is clearly 
unacceptable. This again shows that only T> + can have a solution. In particular, the 
supersymmetric bound state are trustworthy even though the quantum mechanics 
itself would be corrected, at small r, by non-Abelian nature of such objects. 

Therefore, the index problem of the above operator is on par with that of zero 
mode problems around non-Abelian monopoles; the Callias index theorem [42, 29, 
30] should apply. We thus anticipate that the number of zero energy bound states 
is additive; when the core state is composed of many centers of charges 7 C)J 4 with 
(lc,A,lh)fco > 0, the number of the bound state of the probe dyon is the naive one, 




(4.22) 



A 



since 7 C = J2a 1c,a- 
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5 Wall-Crossing from Moduli Dynamics 



5.1 Primitive Wall-Crossing: 7 C + 7^ 

So far, we ignored the precise supermultiplet structures; Our approximation allowed 
us to treat the supermultiplet structure of the core state as a separate sector, while we 
extracted only partial sector of the probe dyons which would have been responsible 
for building a half-hypermultiplet. More generally, the probe dyon can come with 
higher spin states, such as iV = 2 vector multiplet or higher, so we may decompose 
the Hilbert space of the combined core-probe system as 

core ® ^probc ® ^moduli dynamics • 

(5.1) 

The reduced Hilbert space denotes part of the free Hilbert space of a BPS particle 
that multiplies the half-hypermultiplet, 

ft = ft reduced ®([l/2]©2[0]) . (5.2) 

When the probe dyon is in the half-hypermultiplet,* 9 %probc Cd w °uld have only one 
state, while in the vector multiplet, it would be the angular momentum 1/2 Hilbert 
space, etc. 

The decomposition (5.1) can be understood easily. The core part of the Hilbert 
space is inert, so can be treated as non-dynamical. Of the probe, the half-hypermultiplet 
part are generated by the universal would-be Goldstino modes which become no 
longer free due to the presence of the core state. Instead they participate in the 
moduli dynamics we constructed and thus belong to "H mo duii dynamics- Note that these 
four modes would become free at r = 00, regaining its nature as Goldstino. The 
remaining part %probe ed accounts for extra degeneracies and spin content of the probe 
supermultiplet, which should represent additional structure on top of the low energy 
dynamics. 

On the other hand, the second helicity trace (1.3), which is the relevant index for 
N = 2 theories, takes value 

fi([j]®([l/2]©2[0])) = (-l)^(2 J + l) (5.3) 

for the irreducible angular momentum multiplet [j], and can also be expressed as 

il (U) = traduced (-l f h . (5.4) 



#9 Recall that usual hypermultiplet forms when the CTP conjugate states are taken into account. 
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The degrees of freedom for the core state does not participate in the dynamics, so we 
have the decomposition 



^ (^corc ® "^probc ® ^moduli dynamics) 

= {l(H corc ) x traced (-1)*» x tr Wmodulidy _(-l) 2J 3 

= n (n corc ) x n (^ probc ) x tr Wmoduli dynamics (-i) 2 ^ . (5.5) 

Combining with the supersymmetric bound state we found above, this reproduces 
the primitive wall-crossing formula of Denef, 

Afi( 7c + 7 /0 = -(-1) 2 K^>I 2|( 7c , 7 ,)| n( lc )n( lh ) . (5.6) 
5.2 Semi-Primitive Wall Crossing: 7c + n^h 

The semi-primitive wall-crossing formula of Denef and Moore conjectures how many 
BPS states of charge 7c + w-fh appears across a MSW, for positive integer n In order 
to compute the degeneracies of such states we must consider n number of 7^ charges 
in the core state background of 7 C . The Lagrangian would be 

n 

C = J2^) + ^hh , (5.7) 

i=i 

where denotes the one-particle Lagrangian for i-th probe dyon, all of which are 
of the identical form. Chh captures the interaction among (identical) probe particles. 

In our approximation, the latter can be ignored as long as the charges are such 
that 

K7c,7fc>l»l<7fc,7 , fc >l • (5-8) 

In particular this is the case if the probe charges are all mutually local, e.g., the 
same or proportional to each other. Then, the latter term Chh represent the second 
order correction to the former's first order form and can be safely ignored. The only 
nontrivial remnant is the matter of statistics, as in any quantum mechanics of many 
identical particles. 

In addition, there is also a logical possibility that one-particle BPS states of non- 
primitive charge k'jh exist. In supergravity, such states are always there, since black 
holes can have any quantized charges. In field theory setting, the situation is a little 
unclear. In five dimensions, multi-instanton bound state do exist in the maximally 
supersymmetric Yang-Mills theory as quantum one-particle states. However, they 
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are tied to the UV completion of this theory which is the mysterious (2, 0) theories. 
In the more familiar four-dimensional Yang-Mills setting, we are yet to see such an 
example. Nevertheless, we will include the possibility that the probe dyon of our 
moduli dynamics is non-primitive. Then, counting the degeneracy of the bound 
states 7 C + n^h is basically identical to partition of wjh into identical halo particles 
of rry/j = (X]j m i^)7/i with some cares on the statistics of each dyon of charge k^h- 
If it turns out that such non-primitive states do not exist,* 10 we may simply set 
Vt(k lh ) = for k > 2. 

The question of statistics lead us to consider the intrinsic spin of the individual 
probe particle in the moduli dynamics. While the quantum mechanics by itself won't 
tell us about statistics of the particle, we can invoke the usual spin-statistics relation 
and instead ask about the spin. Recall that the canonical commutators, 

{$ m ,ij, n } = 8 mn , (5.9) 

implies that the spatial rotation generators of SU(2)l acting on the wavefunction are 
~\ r , I? 6 ] - \ e abc {r, ^} = \ e abc ( ° £ ) . (5.10) 

This shows that the 4-component wavefunction, consists of a single spin doublet 
V in the lower half and a pair of spin singlet states combined into the upper half part, 
U. Recall that the bound states can appear only in the V sector; the supercharge 
Qa is effectively positively definite on U as we saw in section 4.3. Therefore the BPS 
bound state of a (half-)hypermultiplet probe and the core always involve of a spin 
1/2 wavefunction. 

More generally, the probe might be in a bigger multiplet, where %probe Cd is a l so 
part of the data that enters the probe dynamics although we simply factored it out. 
Taking into account the latter, we can see that the probe particle can be seen as a 
particle of spin content in the moduli quantum mechanics 

Krob C e ed ®([l/2]e2[0]), (5.11) 

but the BPS bound state appears only in the sector "Hprobc Cd ® [1/2]- For example, if 
^probc Cd = total spin of the probe dyon that is involved in the bound state 

formation is S ± 1/2. Therefore, as far as supersymmetric bound state formation 

# 10 The result of the previous section is suggestive in this regard. The bound states exist only if 
the Schwinger product of the two constituent charges are nonzero. Even if we take into account the 
finite core mass, we expect that a single-particle bound state of type £7 + 7 probably does not exist, 
which in induction suggests absence of the state of charge kj for k > 2 altogether. An interesting 
question is how this feature is modified in the realm of supergravity, where black holes of large 
non-primitive charges appears. 
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goes, that the probe dyon can be treated as if it is Boson or Fermion for 2S odd or 
even, respectively. 

Such assignment of statistics is precisely what we expect on the field theory 
ground: Note that 5 = correspond to the hypermultiplet while S — 1/2 to the 
vector multiplet. When one construct BPS dyons in the weakly coupled theory, the 
simplest method is to excite massive electrically charged and L 2 -normalizable modes 
around magnetic soliton [8]. When the charged field is in the hypermultiplet, the rel- 
evant excitations arise all from the Dirac field and the Fermi statistics rule when we 
try to construct the dyons. For a vector multiplet, additional modes arise both from 
the vector field, so the Bosonic statistics become dominant. This naive construction 
works verbatim for N = 4 Yang-Mills theories, while for N = 2 only slightly mod- 
ified (i.e., degeneracy shift by unit) as seen from more rigorous index computation 
[14, 17]. When we phrase the N = 2 result in terms of vector multiplet contributions 
vs. hypermultiplet contributions, we see the above statistics assignment emerging. 

Interestingly, this statistics is correlated with the sign of index Q of the probe 
dyon since 



[S](8)([l/2]©2[0]) = (-1) 2S (2S + 1) . (5.12) 



Thus, in the context of our probe moduli dynamics, probe dyons with positive Q 
should behave as Fermions, while probe dyons with negative Q should behave as 
Bosons. More generally, %probc Cd can be a direct sum of more than one spin sectors. 
We write 

= ®a[Sa] = R+ © R- , (5.13) 

with R± denoting the decomposition according to the sign (— l) 2Sa . Thus, 

^probe = dimR + — dimR_ . (5-14) 

For the purpose of the moduli quantum mechanics here, then, we effectively have 
dimR + Fermions and dimR_ Bosons of the same probe charge. 

Once this statistics issue is cleared, one can construct the generating function for 
the index Q(j c + njh) as follows 

oo 

J2^(lc + n lh )q n = Q( lc )-Tr[(-l) 2J3 q N ] . (5.15) 

n=0 

We used here the notation Tr to emphasize that it is performed also over the dyons 
of various charges kjh as well as over the individual Fock space with the number 
operator N that counts the multiple probe dyons of the same charge. Let us split the 
number operator N = V\ „. 3 „- 3 kNP. 3 ., + Y\ . 3 . 3 kN^ . 3 ., with N B for bosons 
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and N F for fermions. Here |j| xt | < |(7 C , k^ h } \ — \ and < S a . The relevant trace 
then becomes 

h N 



t[(-i) 2 - 



N 



N. 



B/F 

fc ' J ext' J CT 



which can be summed explicitly as 



(oo 

= n [ i ~ (-i) 2fc|<7c ' 7,i>| (? A: 

= n 



2k\( lc , lh )\ q k 



N 



- dim(j ext )-(dim(R + )-dim(R_)) 

2\{ lc ,k lh )\n(k lh ) 



2k\( lc , lh )\ q k 



N 



(5.16) 



It shows that the generating function is 

£ ^ + ^)g n = n(7c) II [ x - (-i) 2fe<7h ' 7c Y 



2fc|(7 h , 7c )|n(fc7h) 



(5.17) 



n=0 



fc=l 



This is precisely the semi-crossing wall-crossing formula conjectured by Denef and 
Moore [18], provided that the one-particle states of charge ^c + nih are absent on the 
other side of the wall. Note that the latter assumption is guaranteed by our moduli 
dynamics. Thus, by staying near the walls of marginal stability and adjusting the 
probe dyon to be much lighter than the core, we have derived the semi-primitive 
wall-crossing formulae from the first principle. 



6 Conclusion and Discussion 

We have derived a M = 4 supersymmetry low energy dynamics that govern probe 
dyons interacting with relatively heavy core states, in the long distance approxima- 
tion. The proximity of the Coulomb vacuum to the marginal stability wall acts as 
a crucial control parameter that allows this non-relativistic quantum mechanical de- 
scription, and we were able to reproduce the conjectured primitive and semi-primitive 
wall-crossing formulae for Seiberg-Witten theory dyons. 
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An important technological step here was to incorporate the potential energy of 
the probe particles, due to the core state, into the supersymmetric quantum mechan- 
ics. Because the latter comes with different bosonic and fermionic degrees of freedom, 
a nonconventional form of the supersymmetric low energy theory emerged, but in a 
manner consistent with the BPS structure of the underlying N = 2 field theory in 
question. 

As we mentioned early on, our approximation scheme was inspired by the notion 
of framed BPS state in presence of a line operator. See Appendix C for a short review 
on line operator in relation to the wall-crossing. In a sense the line operator provides 
a setting where our computation would become an exact description and can aid 
evaluation of the line operator expectation values. The vacuum expectation of line 
operator is in effect a (— 1) F weighted trace over the Hilbert space with a particular 
charge object T inserted as an external object, 

H={Ql,Qc}, (6.1) 

where denote the supercharges preserved by the line operator. It was conjectured 
that this observable can be expanded into 

(L r ), h = J2 n ( T + ^)X, h , (6.2) 

where A^'s are the Darboux coordinates of [44]. The semi-classical analysis on the 
conjectured form of (Lr) would be interesting and illuminating as in Ref. [43]. As 
noted by Gaiotto et.al [44, 45], this asserts the much needed continuity property of 
A"s over the vacuum moduli space that plays a central role justifying KS formalism 
in the context of iV = 2 Seiberg-Witten theory. Our low energy quantum mechanics 
is consistent with this claim since 

n(r + 7fc ,c)A^(c,i2) 

= e- 2nKcK ~ lz ^ ] tr r+gh [(-l) F e- 2 * RH ™^ ~ ieQ a(Q)} x (•••) , (6.3) 

where the first two terms follows from discussions in section 2, while <r(Q) denotes the 
quadratic refinement, as argued in Ref. [20]. The trace is over the quantum mechanical 
Hilbert space for the charge T + jh, while the parenthesis denotes subleading loop 
contribution in the given charge sector. 

An important generalization of our analysis is to study the wall-crossing phenom- 
ena in the N = 2 supergravity. In fact, the formalism we developed is more natural for 
the supergravity system, since the horizon provides natural cut-off at short distance 
and renders the Abelian description of the core state exact. That is, one can hide 



(Li 



Tr 



[-Ife 



-2ttRH 
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the any potential subtlety associated with the Coulombic centers behind the horizon. 
Quantum mechanical description of more than one extremally charged black hole has 
been studied previously, but only in the context of same charge black holes, which is 
a particular limit of our dynamics without potential terms. We are poised to consider 
many black holes with mutually non-local and interacting center, and elevate Denef 's 
old discussion black hole halos to fully quantum level. 

In both field theory and the supergravity version of such a low energy quantum 
mechanics, there is a simpler way to count bound states. As long as the true moduli 
space defined by /C = is compact, the relevant supercharge would be Fredholm, 
and one could compute the index by concentrating on the true moduli space defined 
by K = 0. The quantum mechanics then would reduce to a supersymmetric Landau 
level problem on a curved 2n dimensional manifold, and can be presumably counted 
by computing the volume of this true moduli space. A similar idea has been recently 
used in [46, 47], but our approach provides a rigorous derivation of such a method 
and thus the precise state counting. Details of this computation will be presented 
elsewhere. 

Finally, though we have focused on the moduli space dynamics of framed BPS 
particles in D = 4 N = 2 supersymmetric gauge theories, our analysis can be poten- 
tially applied to study the wall-crossing phenomena of any supersymmetric theories in 
presence of higher dimensional external objects. One potential application is a study 
of the wall-crossing formulae of the four- dimensional gauge theories in presence of a 
surface operator, which has been conjectured in [48] as a hybrid of 2D Ceccoti-Vafa 
WCF [3] and 4D Kontsevich-Soibelman WCF [19]. Our analysis also would be use- 
ful to study the wall-crossing formulae of two-dimensional M = (2, 2) massive CP™ 
models in relation to that of four-dimensional M = 2 SQCD [49, 50, 51, 52]. 
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Appendix 



A BPS Equation for the Semiclassical Core 



This appendix reviews the BPS equation, of Seiberg-Witten low energy theory, for 
long-range Abelian fields for any given core charges. One can easily read off N = 2 
SUSY variation rules in four dimensions from N — 1 SUSY variation rules in six 
dimensions 

1 

2" 

where A and e are six-dimensional chiral spinors, 



5\a — -Fmn^ m ' 



(A.l) 



r 012345 A A = \a , r 012345 e A = e A . (A.2) 

Here A = 1,2 are the R-symmetry indices. Let us decompose the six-dimensional 
gamma matrices T M as 

' 

a» 



L 2 , 



r 4 = 7c 



r 5 = 7c 



Y 



(A.3) 



where 



7 



0123 



. In the above representation, the gaugino \a can be decomposed 



into Xa = X a A © A^. As usual, a, a denote the 4-D chiral/ ant i- chiral spinor indices. 
One can then rewrite (A.l) as 



5\ aA = ^^a^A + ^ ' ccc^A D »<t> , 



(f) = A 4 + iA 5 . (A.4) 
With Z c = \Z C \ (, the core state configuration should satisfy the following relation 

(A.5) 



or equivalently 
that is, 



(Q A + iC 1 Q A a°)E A ,\B 
-if e B -{B + iE- iC _1 W) - C^Bd t (f) = , 
f - <" X V0 = , d t (f) = . 



(A.6) 



(A.7) 



One quick way to show that ( represents the phase factor of Z c is to look at the 
energy for the configuration (A.7), say, for rank one example: performing the usual 
trick of completing the square with (A.7) in mind, one obtain 



— / d 3 x Imr 

8tt 



B 2 + E 2 + |V0p 



= Re 



(A.8) 



with Z c = P0d(oo) + Q<f){oc). This shows that C _: % = \Z C 
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B More on M = 4 Quantum Mechanics 



Here we present more on M = 4 Lagrangian with conformal -R 3 target manifold. Here, 
we first derive the massless case with curved background and then add potential 
terms, which provides an alternate path to (3.28). Then, we spend some time on 
supercharge operators and quantum Hamiltonian. 



B.l Massless and curved 

First of all, we wish to fill the gap between sections 3.1 and 3.3 with a derivation of 
massless M = 4 theory onto conformally flat R 3 , which turned out to be regarded as a 
special case of theories in Ref . [22] . In next subsection, we demonstrate that how the 
massive Lagrangian of section 3.3 emerges by combining the result of section 3.1 with 
this massless case. Based on the educated guess and group theoretical consideration, 
one possible candidate for Af = 4 SUSY transformation rules are following 

5x a = irf mn e m r , 6ij m = V a mn e n x a + ae^J^dJ^" , (B.l) 

where a will be determined. Here i]^ nn denotes the 't Hooft tensor with the convention 
^i 3 2 = 4a = +1- 

To start, consider a standard kinetic term for flat target manifold, 

£(o) = I/iV + l -f^ m , (B.2) 
whose variation under the M = 4 SUSY transformations is 
5(±fx a x a ) = t - V l n dJe m rx a x a + tf V a mn e m rx a , 

£(^"v m ) = - ^;a a /e ? frr - tfv a mn e m rx a - % -rf mn d h fx a x\ m r 

+ »a/fW + ^CJ-^J^iVTeY • (B.3) 
• One can reorganize the velocity-square terms in (B.3) into 



2 



-d b fe m r V b mn x a - V a mn x b x a = -e eab e ecd x a x c d b fv d mn ^ 



i 



2 



+" ^VpmVnp£eabX a X C d b fe m ^ r 



= - -e eab ■ x a d b f V e np srr - -a^r'djd.frji^r^r 

+ \ar l d a fd a fe mnpq i> m rV^ q ■ (B.4) 
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The first term in the last equality of (B.4) implies that we have to add the 
following term 

S{ + \e ahc x a d h M m J, m r) = + l -e abc x a d b f V c m J^ m r 

~ ^a b cV a pq V c mn d b fe p ^ m r - \e abc x a d b dJ V c m X q e p ^ m r ■ 

(B.5) 



Using the identities of 't Hooft tensor 



e abcVmnVpq ^rnpVnq ^npVmq ^nqVmp ^rnqVnp i 

viqV'mn + V^Vnq + rfprjfqm + VpsVrs^mnr = , (B.6) 

one can massage the second and third terms in (B.5) into followings: 

-\e abcV ; qV c mn d b fe^ m r = X -r? mn d a fe m r ■ - ^«/rf • <*V , 

(B.7) 



and 



= + -^e mnP qd 2 w m rvw q 

- ^x a d a d c f V c pl e^y m re q mni ■ (B.8) 

In summary, one can show that 

5{ + l -e abc x a d b w m j, m r) = + \t abc x a d b M m M m r + l -^ mn dje m r ■ 



(B.9) 



Here one can determine, from the fourth term in second equality of (B.3) and 
third term in (B.9), the value of the coefficient a by 

a = + l - (B.10) 
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Collecting all the results so far, one can have 

5(±fx a x a + l -f^ m + % -e ahc x a d h frf m J, m r) 



B.ll) 



At the end of the day, this gives the massless M = 4 non-linear sigma model 
therefore takes the following form 

£<°> =\fx a x a + % -f1> m i> m + \e ahc x a d h jrt m J, m r 

- ^d 2 je mnpq ^ m rrr + y-\d a ffe mnpq ^rrr , (B.12) 

where the covariant derivative for fermions is defined as 

V t ^ m = i, m + ^e abc x a d b hgf V c mn r . (B.13) 
The above massless Lagrangian is invariant under the M = 4 SUSY transformation 

5x a = i V a mn e m r , 5^ m = r) a mn e n x a + '-e^f-'dj^r ■ (B.14) 
This is the curved space version of (3.2). 



B.2 Massive and curved 

Now we wish to add potential terms to this by twisting the supersymmetry transfor- 
mation rules. From discussion of section 3.2, it is clear that the right thing to do, at 
least in the context of M = 1 supersymmetry, is to shift the fermion transformation 
rule as 

5x a = t V a mn e m r , # m = V a mn e n x a + e m - f (k + -^dJVV^ , (B.15) 

since the last piece multiplying e m is nothing but the on-shell value of the auxiliary 
field b. The corresponding Lagrangian from (3.28) 



c = i/ 



x a x a + # m V^ m - jf mnpq {V 2 f - {d a \ogf) 2 }^rw} (B.16) 



-W a x a + id b W c ^ c + tf^daif-^VrX - ^e ab JCf- 1 d a M b r - ^.K 2 
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is indeed consistent with the above massless one in (B.12). 

To show that this Lagrangian is invariant under this transformation, we split it 
into three parts, C = £ (0) + + £ (2) , as 

£ (0) = \fx a x a + l -f^ m + Z -e abc x a d b f V c mn ^ m r - W a x a 



£ (2) = ~\fK\ (B.17) 

where we introduced K = / _1 /C. £ ( - ^ is already invariant under (B.14), so we have 
only K- dependence pieces in 5£(°\ which is 

1 



12 



5'/ - ^r 1 ^/) 2 ] . (B.18) 



After some tedious computation, we find 

+if 1/2 d a U 1/2 K)x a e m ^ m + k a6c i; a / 1/2 9 6 (/ 1/2 X)77^e m Vn 
+6( 1 -fK*) , (B.19) 

which, combined with 5C^°\ give us 

8(£<® + = 5(^/tf 2 ) = -5£ (2) , (B.20) 

as required. 



B.3 Supercharges and Hamiltonian 

Using the M = 4 supersymmetric variation rules (3.30) , the Nother charges of the 
M = 4 supersymmetry therefore become 

Qm = -V a mn r(Pa + W a ) + ^Lf^a/f + ^«W I W 1] + ^ . (B.21) 
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For completeness, let us check whether the above supercharges give the correct 

supersymmetric transformation rules for bosons and fermions. One can read off from 
the Lagrangian (3.28) the canonical quantization 

[x a ,p b ] = i5 a b , {*p m ,r} = r l 5 mn , [p a ^ m ] = y-'daM™ . (B.22) 
One can show 

{-v a mp r(Pa + w a ),r} 
= -v a mn r\p a + w ) - faj-'djrr 

= ~V a mn i a + y- l dafri a n P i> [m r ] - lr 2 d a fv a mn , (B.23) 
where we used for the second equality the definition of momentum operator p a 

Pa + W a = fx a + '-e^Af^^r . (B.24) 
One can also show that 

= ^-r'dj^rr + \r x d a w nv ^ m + \r 2 d a w mn , (B.25) 

where we used an identity of 't Hooft tensor 

^abcVLnVpq = Vmp S ng ~ Vnp S mq + Vnq S mp ~ Vmq S np ■ (B.26) 

It implies that 

{Qm, a,} = -r] a mn x a + <w _1 [k + -j-'dj^rr^ , (B.27) 

while the action of supercharges on the bosons follows immediately, 

[Q m x a ] = <„f . (B.28) 

These are precisely the supersymmetry transformation rules in (3.30). 

Finally, we wish to determine the quantum form of the Hamiltonian using 

Ql = H . 
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Let us first write 

Q 4 = V°(p + W) a + A(/C + Z) , 

where 

Z= l -d a fr^ + \e abc d a ffy c . 

Using {Q 4 , A} = (/C + Z)/f and {Q/t,^"} = i a = / _1 vr a , with the supercovariant 
momentum operator 

= (P + W)« + r o , r a = 1 d b f^ b r ] - \ eabcdbf^ , 

we find 

{Q 4 ,Q 4 } = {Q 4 ,^ a (p + W) a + A(/C + Z)} 

= j7T a ( P + W) a + j(IC + Z) 2 

-riQ*, (p + W) a ] - A[Q 4 , £ + Z] . (B.29) 

Let us separate out terms involving either W or /C from the last two terms. Using 
d/C = *g?W, we find 

{Q 4 ,Q 4 } = jn a (p + W) a + j()C + Z) 2 -2id a lCr\-iea b cdaJCiJ b iJ c 
+ (ri(p + W)„, + A[Z, /]) (p + W) 6 
+ (^[(p + W) a ,A] + A[Z,A])Z 

+2^ a A[(p + W) a ,Z] . (B.30) 
By explicit computation one can see that 

{r[(p + w) a ^ b ] + \[z^ b }) = jT b 

(^ a [(p + >V) a ,A]+A[Z,A]) = 0. (B.31) 
Since [(p + W) a , r o ] = upon the summation over a, 

* "( p + W ) a + Ir b ( p + >v) b = 1^- 1 r a r a . (B.32) 
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Finally expanding (/C + Z) 2 out, we complete the potential terms associated with /C 
from K? + 2/CZ, but have a leftover piece Z 2 . So combining them all, we have 

+jZ 2 - - r h r b + 2rx[(p + w) a , z] (B.33) 

The last line can be organized in terms of the curvature of the fermion bundle, 

[D a , D b ] = F abmn ij m r, D a = d a + iT a , (B.34) 
and has the explicit form, 

-\F ahmn r^ m r = ^{2(d 2 f)-f-\df) 2 )e mnk ^ m r^ 1 

-\r\d 2 f)+ l -f-\df) 2 , (B.35) 
Thus, the Hamiltonian H = {Q^,Q4}/2 is 

H = jjn a n a - \F ahmn r^ m r + ^jK 2 - % - ^J^d^r^WV (B.36) 

Although SU(2)r is not manifest in the curvature piece, it is actually SU(2)r in- 
variant as can be seen from (B.35). This coincides with the classical Hamiltonian up 
to normal ordering; the curvature pieces generate extra terms because quantum ^'s 
obey not the Grassman algebra but the Clifford algebra. 

Note that the kinetic term is slightly unconventional in its choice of normal order- 
ing. Because of this, the inner product in the Hilbert space of this quantum mechanics 
should be defined as 




(B.37) 



More usual choice of kinetic term/inner product is related to our convention by rescal- 
ing of the wavefunction by a factor of / 1//4 . 



C Review of KS Invariant and Line Operator 

The idea of the framed BPS state originally arises in study of four-dimensional N = 2 
supersymmetric theories in presence of an external particle of charge T, called line 
operator L T . The line operator can be characterized by the phase factor ( of its 
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central charge Z r . Compactifying the theory on a circle, it has been conjectured 
in [20] that the vacuum expectation value of L r can be expanded in terms of the 
Darboux coordinates X 1 with integer coefficients 

(L r ) = ^fi(r + 7 )A' 7 , (C.f) 

7 

which provides us a direct physical interpretation of Darboux coordinates. Each 
integer coefficient Q(F + 7) here represents the supersymmetric index of a framed 
BPS state of charge 7 bounded to Lp. The Darboux coordinates are very useful to 
compute the hyperKahler metric on the Coulomb branch of four-dimensional theories 
on a circle. 

The expectation value of the line operator depends on both ( and the Coulomb 
branch parameter a in four-dimensional theories. Due to the fact that the physical 
observable (Lr) should not have any discontinuities as ( and a change, important 
consequences of (C.l) are that one can understand how the Kontsevich-Soibelman 
invariant naturally arises, and that provides the origin of the thermodynamic Bethe 
ansatz equation the Darboux coordinates should satisfy. 

Let us now review in this section the central importance of semi-primitive wall- 
crossing formula to derive the Kontsevich-Soibelman BPS invariant in the context of 
line operators. For more details, it is referred to [20] . 

As discussed in the main context, the Witten index Q(T + 7, () can jump once 
the phase of central charge for a certain probe(halo) particle of 7^ is parallel to that 
of the external particle of T denoted by arg(£). That is, when ( moves across the 
so-called BPS ray hi = {( \ Z h /( e the index could have discontinuity. One 

advantage on computation of the index jump in presence of line operator is that the 
wall-crossing phenomena is essentially restricted to the semi-primitive ones. 

Let us now consider the vacuum expectation value of the line operator conjectured 
as in (C.l) 

(L r ) = ^fi(r + 7 )* 7 , 

7 

where X 1 satisfy a multiplication rule below 

X^X 12 = (-l) 2 <^>* 71+72 . (C.2) 

Let us then increase the phase parameter arg(£) so that it moves across the BPS ray 
k- 

• Look at the relation (2.22). If (7 C , 7^) > 0, we have a stable bound state between 
core and halo particles before ( cross the BPS ray k- Then, one can reorganize 
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(C.2) before across the ray into the following form 

(L r )_ = Y, X ^-Y.^ + n ^- l ) 2nhcnh)X l ' 

2n(7c7hM™7h) 



7c 



n=0 



n=l 



(C.3) 



Note that we used the semi-primitive wall-crossing formula (5.17) for the last 
equality. Since we loose the Fock space of halo particles after across the ray 4, 
one can say that 



(Li 



(C.4) 



One can therefore conclude that, since (L r ) should be continuous across the 
ray, X lc is required to jump across the wall by the amount 



X lc ->• X Jc 



nf 1 -^ 



n=l 



2n(7h.7c)f^(«7h) 



= n /c Sr ) (^7 C ) ? (c.5) 



where 



^7h('^7c) ~~ ^7 



1 - X 



lh 



n=l 



2<7fc,7c> 



(C.6) 



It is noteworthy here that this is the desired discontinuity how the Darboux 
coordinate X^ jumps across the BPS ray lh- 

Let us now in turn consider the converse, i.e., (7 C , 7h) < 0. According to (2.22), 
there is no stable bound state between the core and halo particle before the ( 
across the BPS ray l h . Then, one can rewrite (C.2) before across the ray as 



(It)- - ^ & lc X lc . 



(C.7) 



7c 



Since we gain the Fock space of halo particles after across the ray l h , one can 
say that 



(L r )+ = J>(7c)* 7c J] [l - * 



7c 



n=l 



-2n(7 c ,7 h >n(n7 h ) 



(C.8) 



X lc is again required to jump across the wall by the same amount 



*7 C *7 C n [ i - 

n=l 



X, 



lh 



2n(7 h ,7 c )n(n7, l ) 



II^SrH^c), (C9) 



n=l 



which is the same to (C.6). 
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Let us now consider two chambers of A^couiomb x C*, the Coulomb branch and £- 
plane, separated by walls of marginal stability. The physical observable (L r ) should 
not depend on choice of a path connecting those two chambers. The different paths 
however in general cross different set of walls of marginal stability. One can therefore 
conclude, from the fact that there are infinitely many possible line operators, that a 
path-ordered product of transformations below 

I =UU IC nt h) (CIO) 

lh n 

defines an invariant over the Coulomb branch A^couiomb- is indeed the so-called 
Kontsevich-Soibelman invariant. 
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